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ABSTRACT 

We study the runaway mass loss process of major eruptions of luminous blue 
variables (LBVs) stars, such as the 1837-1856 Great Eruption of rj Carinae. We 
follow the evolution of a massive star with a spherical stellar evolution numerical 
code. After the star exhausted most of the hydrogen in the core and had de- 
veloped a large envelope, we remove mass at a rate of 1M© yr _1 from the outer 
envelope for 20 years. We find that after removing a small amount of mass at a 
high rate, the star contracts and releases a huge amount of gravitational energy. 
We suggest that this energy can sustain the high mass loss rate. The triggering 
of this runaway mass loss process might be a close stellar companion or internal 
structural changes. We show that a strong magnetic field region can be built in 
the radiative zone above the convective core of the evolved massive star. When 
this magnetic energy is released it might trigger a fast removal of mass, and by 
that trigger an eruption. Namely, LBV major eruptions might be triggered by 
magnetic activity cycles. The prediction is that LBV stars that experience major 
eruptions should be found to have a close companion and/or have signatures of 
strong magnetic activity during or after the eruption. 



1. INTRODUCTION 

Luminous Blue Variables (LBVs) are massive hot luminous stars. They posses very 
strong winds that exhibit irregular variabilities on time scales ranging from days to years. 
On top of these variations, LBVs experience extreme mass loss rate episodes (e.g., Smith 
& Owocki 2006; Owocki & van Marie 2009, and references therein), e.g., the 19th century 
eruptions of rj Carinae (Humphreys et al. 1999), where a mass of ~ 10 — 2OM was lost 
(Smith et al. 2003b; Smith 2006; Smith & Owocki 2006; Smith & Ferland 2007). These 
eruptions cannot be accounted for by the regular stellar luminosity, and they require some 
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extra energy source, e.g., internal structural change in the star, that might even increase 
the stellar luminosity above the Eddington limit (Owocki & van Marie 2009). However, 
part, or even all, of the increase in the luminosity of rj Car in the 1837-1856 Great Eruption 
could have come from gravitational energy of the mass accreted by the secondary star (Soker 
2007). The accretion of mass onto the secondary star can explain also the kinetic energy 
of the Homunculus (Soker 2007); the Homunculus is the bipolar nebula of r] Car that was 
formed in the Great Eruption (Davidson & Humphreys 1997). 

The influence of radiation on the mass loss process of stars near their Eddington lu- 
minosity limit in relation to LBV eruptions is discussed by van Marie et al. (2008, 2009) 
and Owocki & van Marie (2009). In particular, they discuss how the extended porosity 
formalism (Shaviv 1998, 2000) can account for the Great Eruption of rj Car (Owocki et al. 
2004; van Marie et al. 2008, 2009). In the present paper we do not deal with the interaction 
of radiation with matter. We rather limit ourself to discuss the possible instability process 
that can lead to the release of a huge amount of energy by internal structural change. 

One of the significant differences between our approach and most other studies of the 
Great Eruption of rj Car concerns the energy of the Homunculus. While most studies (e.g., 
Smith 2006) attribute the entire energy source to the primary enhanced luminosity, we take 
the view that most of the kinetic energy of the Homunculus results from two opposite jets 
that were blown by the companion during the Great Eruption. The companion blew the 
jets as it accreted mass from the primary dense wind via an accretion disk (Soker 2007). 
Therefore, we do not deal with the energy of the homunculus, but only with the energy that 
is required to unbind a mass of ~ 10 — 2OM during the Great Eruption of r] Car, and similar 
LBV eruptions. The same jets can account for fast ejecta that were blown from the rj Car 
binary system (Smith & Morse 2004). 

The mass of ~ 10 — 2OM (Smith et al. 2003b; Smith & Ferland 2007) that was ejected 
in the Great Eruption resides in the outer part of the radiative outer region of LBV stars. 
In section 2 we build a stellar model that has a similar structure to that of rj Car before the 
Great Eruption, and discuss some of its properties. Soker (2007) already speculated that 
the Great Eruption of i] Car was triggered by disturbances in the outer boundary of the 
inner convective region, most likely by magnetic activity, that expelled the outer radiative 
zone. Soker (2007) further mentioned that one way to form an extended envelope is by the 
contraction of the inner layers. In section 3 we go one step further and show that indeed, the 
removal of the outer region of the star causes the star to shrink and release a huge amount 
of gravitational energy. Earlier suggestions for the cause of LBV instability was summarized 
by Humphreys & Davidson (1994). Some of them, e.g., the geyser model of Maeder (1992), 
cannot work for a hot star like rj Car. As we show, our model works for blue stars. In section 
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4 we discuss in more detail the possibility that the initial mass removal in LBV outbursts is 
triggered by magnetic activity. We summarize in section 5. 

2. STELLAR STRUCTURE 

We evolve a spherical stellar model with the same evolutionary code that was used 
by us in previous studies over the years (for detail see Soker & Harpaz 1999). We start 
at t — with a zero- age main sequence star of mass Mo = 190M Q . Mass loss is not a 
major part of our study as we are interested in the stellar model toward the end of the main 
sequence. We simply set the mass loss rate to be M = 2 x 10~ 5 M Q yr _1 (for more detail 
on the evolutionary track of massive stars the reader can consult, e.g., Meynet & Maeder 
2003, 2005). The mass, luminosity, and effective temperature, at 4 evolutionary points 
are [M(M ), L(1O 6 L ), T e (10 4 A)] = (190,3,5.7), (160,3,5.1), (150, 3, 4.4), and (139,3,1.6). 
The luminosity does not evolve much, but as the hydrogen in the core is close to exhaustion 
the envelope swells and the effective temperature decreases (see also Smith & Conti 2008). 
In Fig. [I] we show the stellar structure at t — and at t = 2.55 Myr. 

Most relevant to us is the entropy profile. The regions where the entropy profile is 
flat (actually decreasing very slowly) are convective regions. At early times the star is 
almost completely convective. At later times the entropy is flat in the inner ~ 80M Q . The 
outer regions are mainly radiative. Above the inner convective region the entropy increases 
substantially with mass (and radius). Then, in the outer ~ 15 — 2OM the profile becomes 
shallow, and a second convective region exists there. At late times most of the volume of 
the envelope is an outer extended region with very low density (~ 10~ 7 - 10" 6 e cm ) that 
contains a relatively small amount of mass (< 1M & ). 

The evolutionary numerical code calculates the entropy S e using the full equation of 
state. To further elaborate on the entropy behavior to be used later, we examine the quantity 
Sj = Pp 7ad . As evident from Fig. [2J in the massive radiative region above the convective 
core a value of 7 a d = 1.33 accurately describes the rapid entropy rise. In Fig. [2] we plot the 
logarithm of the pressure and of the density, the mass, the accurate entropy calculate by the 
stellar code S e , and of S 1 for 7 a( j = 4/3 (units are given in the caption), as function of stellar 
radius for the second model shown in Fig. [TJ We note again the rapid rise in the entropy 
from the core, and then the flattening in the outer ~ 15 — 20M Q , where a second convective 
region resides in the region 22 < r < 28R & . The very extended outer region is not shown, 
as it contains a small amount of mass < 1M . 

The important property to take from the graphs is that the mass expelled in eruptions 
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Fig. 1. — Density (in units of 0.1 g cm~ 3 ), temperature (in units of 10 6 K), entropy S e (in 
relative units), and radius (in solar radii), as function of mass for two evolutionary times: 
t = on the upper panel and t = 2.55 x 10 6 yr in the lower panel. The two models have 
[M(M ),L(lO 6 L ),T e (lO 4 ^)] = (190,3,5.7), and (139,3,1.6), respectively. The thick lines 
on the horizontal axis mark the convective regions. In the lower panel the photospheric 
radius is 22Oi? , and it is outside the graph. 
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Fig. 2. — Logarithm of the density, logp(g cm -3 ) + 10, of the pressure, in c.g.s units, mass 
in units of 1OM , scaled accurate entropy from the stellar numerical code S e , and 5 7 = 
3.0 x 10~ 16 Pp~ 4 / 3 in c.g.s. units, as function of radius in solar radii. There are two convective 
regions marked by thick lines on the horizontal axis. One is in the core r < 6R Q (M < 
8OM ), and the other at 22 < r < 28i? (127 < M < 132M ). the photosphere is at 
a radius of i? p h = 220R Q , but the region above r = lO5i? contains very little mass, and 
is not shown in this figure. The model is the one drawn in the lower panel of Fig. [Q 
[M(M ),L(lO 6 L ),T e (lO 4 iT)] = (139,3,1.6). 
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of LBVs, such as the Great Eruption of 77 Car, is a high-entropy gas. 

3. THE ERUPTION PHASE 

Stars with a radiative envelope shrink as they loss mass on a time scale shorter than 
the thermal time scale (Webbink 1976; Heisler & Alcock 1986; Maeder 1992). The release 
of gravitational energy by the contracting envelope can lead to an increase in the mass loss 
rate, resulting in a runaway mass loss process. 

To examine the behavior of our model we start with the star at the evolutionary point 
[M{M Q ), L{10 6 L & ),T e {10 4 K)] = (139,3,1.6), and we remove AM burst = 20M o with a con- 
stant mass loss rate of Mb urs t = IMq yr -1 for 20 years. The mass is removed from the outer 
radiative region. This mass loss rate mimics the average mass loss rate during the Great 
Eruption. We then reduce the mass loss rate to M p = 2 x 1O~ 4 M yr _1 , and follow the star 
for another 200 years. In Fig. [3] we plot the radius and luminosity of the star as function of 
time during the eruption. 

The mass loss rate during the 20 years eruption proceeds on a time scale much longer 
than the dynamical time scale, but it is shorter than the thermal time scale. At the beginning 
of the outburst the stellar radius is R = 220-R©. The average thermal time scale of the outer 
region of mass dm is r t h = GMdmj RL, while the mass loss time scale is r m i = dm/ M. Their 
ratio is 

sl =7 (-*l.)( m V-g-vv — y. (i) 

The thermal time scale is substantially longer than the mass loss time scale. As a result of 
this the star losses its thermal equilibrium and rapidly contracts, i.e., on a time scale of few 
years which is much shorter than the thermal time scale. As our model is not fully built to 
take into account evolution on time scales shorter than the thermal time scale, e.g., it does 
not take into account the energy required to remove the mass from the envelope, our results 
of the rapid mass loss episode are not accurate. The code cannot handle properly the removal 
of mass during a time shorter than the thermal time scale, and there are large fluctuations 
in the luminosity and radius. In reality, the star exceeds its Eddington luminosity limit and 
is expected to be unstable (Shaviv 2001). However, as the mass loss time scale is longer than 
the dynamical time scale, our results correctly give the general description. 

It is evident that as mass loss starts on a timescale shorter than the thermal time scale 
the inner regions of the star contracts, and release a huge amount of gravitational energy. In 
our non-dynamical model the released gravitational energy leads to a higher luminosity. In 
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Fig. 3. — The radius and luminosity of the star during the runaway mass loss episode 
(eruption), starting two years before the high mass loss rate and ending two years after. The 
pre-outburst radius in our model is ~ 2OOi? . During the eruption phase, < t e < 20 yr, 
mass was removed from the outer radiative zone at a rate of Mb urst = 1M yr -1 . Because the 
mass loss time scale is shorter than the thermal time scale, the numerical code is unstable, 
and there are large fluctuations. However, the general shrinkage of the star and the release 
of a huge amount of energy is evident. 
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reality, the increased luminosity will result in a higher mass loss rate, even if with a smaller 
terminal speed than in the regular wind, because the luminosity exceeds the Eddington limit 
luminosity by more than an order of magnitude. In particular, the opacity of the expelled 
mass is huge, and the extra luminosity would be absorbed in the wind, and accelerates it. 
Our model does not have the wind component. Therefore, the luminosity given above is not 
the observed luminosity, but rather the energy available to accelerate the expelled mass. 

The total kinetic energy of the gas in the Homunculus is ~ 5 x 1O 49 (M#/2OM ) erg 
where Mjj is the mass in the Homunculus (Smith 2006; Soker 2007). The total gravitational 
energy liberated by the contracting star is ~ 5 x 10 50 erg, which is much larger. However, 
as mentioned, our numerical code does not have in it the energy that is required to remove 
the mass. After losing only 4M the stellar radius has shrunk to ~ 15i? . Losing another 
16M from an average radius of 15i? requires an energy of ~ 5 x 10 50 erg. Namely, most 
of the energy will go to unbind the expelled gas, and accelerate it. As mentioned earlier, in 
the binary model (Soker 2001, 2007) the accreting companion can supply most of the energy 
of the Homunculus, and energy considerations is not of a worry. 

Our model is not fully consistent in another manner. The numerical code does not 
include the triggering process of the high mass loss rate during the eruption. For example, the 
triggering can come from a tidal interaction with the companion which deposited energy to 
the envelope. Tidal enhancement of mass loss rate during the Great Eruption was mentioned 
by Soker (2001), and during present periastron passages by Smith et al. (2003a). Indeed, 
at the onset of the eruption the stellar radius was ~ 1AU, while the periastron distance in 
i] Car is ~ 1.5 AU, implying that a strong tidal interaction took place before and during 
the onset of the eruption. The primary star could have gone through a weak instability 
(pulsation) that increased its radius a little. If this continued while the secondary star was 
at periastron passage, then the secondary could have remove some amount of mass, and 
trigger the outburst. Another possible triggering process is discussed in the next section. 

Although the photospheric radius decreases by more than an order of magnitude during 
the eruption, mass shells do not contract much. The main contraction of the photosphere is 
due to removal of the outer high entropy layers. On average, during the 20 years eruption 
the mass shells at the photosphere (each time a new mass shell) have their radii smaller by a 
factor of ( r < 2 than their radii at the beginning of the eruption. This implies that the ratio 
of rotation velocity to break-up velocity of the mass shells increases by an average factor of 

1/2 

only ( r — 1.5. In any case, even after contraction the mass shells at the photosphere have 
their radii larger than their main sequence value. This suggests that a single star cannot 
posses fast rotation. This discussion further supports the claim that single star models cannot 
explain the bipolar structure of the Homunculus, and an interaction with a companion is 
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required (Soker 2004, 2007). 

After we terminate the high mass removal rate, i.e., the end of the eruption phase, and 
set the mass loss rate to M p = 2 x 1O _4 M yr~ 4 , the star relaxes on a very short time, 
and starts a slow recovery. While during the removal of mass at a high rate the stellar 
layers contract on average and release gravitational energy, on the recovery phase the star 
re-expands. This requires energy, and the luminosity is lower than its pre-eruption value. 
It starts with L = 1.6 x 1O 6 L immediately after the eruption ends, and increases slowly 
to L = 2 x 1O 6 L at 200 years after the eruption. This slow recovery might be related 
to the secular brightness of rj Car. In our model the effective temperature 200 years after 
the eruption is very high, T e ~ 5 x 10 4 K. However, our treatment of the eruption is not 
accurate, and it is possible that the recovery phase is more rapid than what we find. In any 
case, the dense wind of the primary in rj Car blocks the UV radiation from its photosphere 
(Hillier et al. 2001, 2006). 



4. MAGNETIC ACTIVITY 

In §3] we found that an initial rapid mass loss will lead to a runaway process. The 
initial rapid mass loss might be triggered by a tidal interaction with a companion, or from 
internal changes in the star. The proximity of the star to its Eddington luminosity limit 
most probably is also a key property of the star. Here we examine a possible triggering by 
magnetic fields (Soker 2007). 

In the evolved star just prior to eruption, the inner ~ 8OM region is convective (the 
convective core), while most of the envelope is radiative. It is quite possible that a dynamo 
will be operating in the boundary of the convective core and the radiative region above it 
(e.g., Charbonneau & MacGregor 2001; MacGregor & Cassinelli 2003). This boundary is 
located where the entropy starts to rise, at (M, r) = (81M ,6-R ) in the model presented 
in the lower panel of Fig. [T] and in Fig. [2J As can be seen, above this radius the entropy 
rises sharply, until it flatten again in the outer ~ 2OM . The mass of ~ 2OM above the 
steep entropy region, is about equal to the mass expelled in the Great Eruption of r] Car. 
In our model this is not a coincidence. In the sun, magnetic flux tubes that are formed in 
the lower boundary of the convective region buoy outward to the surface. However, in the 
case studied here the magnetic flux tubes need to buoy through the radiative region, where 
entropy increases outward. We now examine this situation. 

We note that a more thorough study is conducted by MacGregor & Cassinelli (2003) 
who considered also stellar rotation and heat transfer. The aim of our simpler treatment is 
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to emphasize the important role that the steep entropy rise in the radiative zone might play 
in forming a strong magnetic field region. 

Consider a magnetic flux tube, or a magnetic flux loop, formed at the bottom of the 
radiative zone r = r , with a magnetic field of intensity B t0 and density p t o, and with an 
initial temperature equals to the ambient (envelope) temperature T t0 = T e (r ) = T e0 at 
this radius. The ambient density is p e , and the magnetic field in the ambient medium is 
assumed to be much weaker than that in the tube (or loop), such that the ambient pressure 
is P e (r) = p e kT e / pyrin, where prriH is the mean mass per particle. Let the cross-section of 
the tube be A(r), with A = A(r Q ), and let its length be L. The flux tube and our stellar 
model just before the eruption are drawn schematically in Fig. HI 

For the evolution of the magnetic field inside the tube we take 

B t = Bto— p = B t o ( — J , (2) 



A \pto 

where the first equality results from magnetic flux conservation, while the value of 5 depends 
on the relative variation of A and L. If the length of the tube does not change much as 
it rises, as might be the case if the magnetic field is strong, then 8—1. If, on the other 
hand, the field is random or the relative expansion of all dimensions is the same, as expected 
when the field is weak and it is dragged by the rising tube, then 5 = 2/3. We will take here 
5 = 2/3 as appropriate for a weak initial magnetic field. 

We also define the initial ratio of thermal to magnetic pressure inside the tube 

P Q — f ^thermal \ _ PW^tO f -Bf \ 1 

The pressure inside the tube includes magnetic and thermal contributions, and can be written 

as 

P - P + P - PtkTt 4- ^ m 
-ft — -fthermal + 1 — ■ (4 ) 

pVTLu 07T 

We assume adiabatic evolution of the tube such that the temperature inside the tube evolves 
as T t oc p7 ad_1 , and use equation (El) to write for the pressure inside the tube 

PtokT to ( Pt y* d + Bl(Pt_\ 2 \ (5) 



/jLTUh \PtoJ »7T \p t0 

The adiabatic index is not constant in the star. However, as evident from Fig. El it is quite 
constant with 7 a( j ~ 4/3 in the radiative zone above the convective core (because radiation 
pressure is very important there). In that zone, 6 < r < 22i? Q (80 < M < 12OM ), the 
entropy steeply rises by a factor of S e (22R Q ) / S e o = 2.3 (Fig. E])- 




Fig. 4. — Schematic drawing of our stellar model just before the eruption, with a flux loop 
embedded in the radiative envelope. There is a core convective region in the sphere r < r = 
6R Q , and a second convective shell in the range 22 < r < 28i? (127 < M < 132M ). 
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Using the definition of (3q we can cast equation (jSJ) into 
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where P e0 = Pto — (1 + P~ x )ptokT t Q/ fimn, as the tube and ambient pressures at the origin 
r are equal. Since we assume that the temperature of the tube is equal to the ambient 
temperature at its origin, we find the expression that connects the density of the tube to 
the ambient pressure at the origin p e0 = (1 + l3~ l )p t Q. As the tube rises, its pressure equals 
to the ambient pressure Pt = P e {r). Using these expressions for p t o and Pt, in equation ([6]) 
gives 
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We now examine the possibility that the flux tube reaches a new equilibrium position, 
namely, a radius where its density is equal to the ambient (envelope) density p t = p e {r). 
This will be possible only if the entropy of the envelope increases fast enough with radius. 



Substituting this equality in equation ([7]), and dividing by [p e (r)/p e0 ] 7ad , we find 
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We take for the envelope entropy S e (r) ~ S^fY) = P e j p^ A , with 7 ad ~ 4/3, and for the 
behavior of the magnetic field 28 ~ 4/3, and derive from equation (jSJ) 
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(9) 



Rising magnetic flux tubes will stop rising when the equality in equation holds. This 
situation is not completely stable. MacGregor & Cassinelli (2003) have showed that when 
heat transfer (by radiative diffusion) and rotation are considered, the flux loops emitted 
in the equatorial plane can attain a stationary equilibrium stable with respect to small 
displacements in radius, but are unstable when perturbed in other directions. However, 
the flux loops move outward on a long time scale (see below). In any case, magnetic field 
energy can attain a large value in the radiative region with the steep entropy rise. The rise in 
entropy in the region 6 < r < 22i? Q (80 < M < 120M Q ) is by a factor of S e {22R Q )/S e0 = 2.3 
(Fig. [2]). From equation (Q we find then, that the magnetic flux tube will come to rest in 
this radiative zone if po ^ 1- Namely, the initial magnetic pressure must be smaller than the 
thermal pressure in the tube. 
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The rapid (on several dynamical time scales, or about a year) release of the magnetic 
energy stored in the radiative region between the two convective regions might trigger a high 
mass loss rate for a short time, that will then sustain itself and develop to the runaway loss 
of > 10M Q from the envelope, as discussed in the previous section. 

If stellar rotation and heat transfer are considered the flux tubes do not reach an equilib- 
rium (MacGregor & Cassinelli 2003). However, the rise of flux loops is determined by heating 
via radiative diffusion. This heating proceeds on a long time scale, > 10 4 yr (MacGregor 
& Cassinelli 2003), which is longer than the expected time between major LBV eruptions 
(Smith & Owocki 2006). Even meridional circulation (MacGregor & Cassinelli 2003) will 
transport the flux loops on a long time scale, much longer than buoyant time scale in a con- 
vective region. In our scenario it is assumed that LBV stars posses magnetic activity cycles, 
and the activity builds itself to a maximum on a time scale shorter than the meridional 
circulation and heat transfer time scales. The magnetic energy that is stored is required to 
remove a large amount of mass in a short time to set the runaway mass loss episode. 

In addition to the amplification of magnetic fields in the boundary of the inner convective 
region and the radiative region above it, the outer convective region can also amplify the 
magnetic field. All these possible dynamo activities are beyond the scope of the present 
paper. Our only goal is to point to the possibility that LBV large eruptions are triggered by 
magnetic activity. 

5. SUMMARY 

Our main goal was to examine the triggering and runaway mass loss process of major 
eruptions of luminous blue variables (LBVs) stars, e.g., the 1837-1856 Great Eruption of r) 
Car. 

We followed the evolution of a massive star. We started with a solar composition and 
with [M(M ),L(lO 6 L ),T e (lO 4 A")] = (190,3,5.7), for the mass, luminosity and effective 
temperature, respectively. We then evolved the star to almost complete depletion of hydro- 
gen in its core (§2), where the star had [M(M ), L(1O 6 L ), T e (WK)} = (139,3,1.6); our 
treatment of the atmosphere is not accurate, and for that the effective temperature can be 
in the range ~ 15, 000 — 20, 000 K. The star has a convective core, then a radiative zone, 
followed by another convective shell (Figs. [Hand [2]). The outer most zone is a very extended 
radiative one, with very small amount of mass. The most important property of the model 
is the development of a steep entropy rise above the convective core. The mass of ~ 2OM 
above the steep entropy region (see lower panel of Fig. [T] and Fig. HJ) is about equal to the 
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mass expelled in the Great Eruption of 77 Car. In our model this is not a coincidence. 

In §3 we studied the response of the star to the eruption phase. We removed 20M Q at 
a rate of 1M Q yr _1 . Numerically, the mass was removed from the outer radiative zone. This 
mass loss time scale is shorter than the thermal time scale of the star (eq. [T|), and before 
removing even 1M , the envelope losses its thermal equilibrium. The envelope contracts, 
and its luminosity increases (Fig. [3]). The increase in luminosity that result from the release 
of gravitational energy is huge. Very likely, this will sustain the high mass loss rate, and 
develop into a runaway mass loss eruption. Our code is not consistent in that we do not 
follow the energy required to remove the mass; we simply remove mass from the envelope. 
Still, we can safely conclude that if a triggering mechanism can remove mass, even a small 
amount of < 0.5M Q , within a time much shorter than the thermal time scale of that mass, 
a runway mass loss process will develop. 

It should be noted that although the photosphere shrinks a lot, the motion of mass shells 
inward is not by a large factor. This implies that the star will not spin-up much during the 
eruption phase. The formation of a bipolar nebulae seems to require a binary companion. 

In §3 we also followed the star for another 200 years after the end of the eruption, with a 
mass loss rate of M = 2 x 10~ 4 to mimic the behavior of rj Car after the Great Eruption. As 
our treatment of the eruption phase is not fully consistent, we also do not have the correct 
values for the post-eruption phase. For example, taking into account the energy required 
to remove the mass from the star would reduce the energy radiated by the star, and would 
make the envelope shrinkage smaller. The recovery phase after the eruption will be shorter. 
Still, we managed to show that the star tries to recover from the eruption. 

The triggering of the rapid mass loss rate can come from a companion or from internal 
structural changes. In the case of r] Car, a strong interaction with the companion took place, 
and this could have triggered the eruption. In §4 we examine the possibility that a magnetic 
activity cycle is the trigger. We showed that a strong magnetic field region can be built in 
the radiative zone above the convective core, as schematically drawn in Fig. HI When this 
energy is released as the magnetic cycle reaches its peak, it can trigger a fast removal of 
mass, and by that trigger the eruption. 

The prediction is that LBV stars that are experiencing major eruptions, should be found 
to have a close companion, and/or have signatures of strong magnetic activity during or after 
the eruption. 

We thank Nathan Smith for useful comments. This research was supported by grants 
from the Israel Science Foundation, and from Asher Space Research Institute at the Technion. 
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